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Nuclear effects for J/^f production in pA collisions are controlled by the 
coherence and color transparency effects. Color transparency onsets when 
the time of formation of the charmonium wave function becomes longer than 
the inter-nucleon spacing. In this energy regime the effective break-up cross 
section for a cc dipole depends on energy and nuclear path length, and agrees 
well with data from fixed target experiments, both in magnitude and energy 
j> ■ dependence. At higher energies of RHIC and LHC coherence in cc pair pro- 

duction leads to charm quark shadowing which is a complement to the high 
i/-) ! twist break up cross section. These two effects explain well with no adjusted 

parameters the magnitude and rapidity dependence of nuclear suppression of 
J/^f observed at RHIC in dAu collisions, while the contribution of leading 
twist gluon shadowing is found to be vanishingly small. A novel mechanism 
of double color filtering for cc dipoles makes nuclei significantly more trans- 
parent in AA compared to pA collisions. This is one of the mechanisms which 
make impossible a model independent " data driven" extrapolation from pA 
to AA. This effect also explains the enhancement of nuclear suppression ob- 
served at forward rapidities in AA collisions at RHIC, what hardly can be 
related to the produced dense medium. J/\& is found to be a clean and sensi- 
tive tool measuring the transport coefficient characterizing the dense matter 
created in AA collisions. RHIC data for px dependence of J/^f production 
in nuclear collisions are well explained with the low value of the transport 
coefficient % < 0.5GeV 2 /fm. 
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1. Time scales and different regimes for J /S& attenuation in nuclei 

The widely used model for high twist nuclear effects is based on unjustified 
assumptions: (i) J/^f (a c-cbar dipole) is always created momentarily inside 
the nucleus; (ii) the produced cc dipole attenuates in the nucleus with a 
break-up cross section a a b s , which is assumed to be universal for all nuclei, 
independent of energy and xf, and is fitted to data. 

Let us start up with examining the latter assumption (ii), assuming for 
the moment that the former one (i) is true. 

1.1. J/'ty formation, color transparency, break-up cross section 

A cc dipole is produced with a small separation r 5c ~ l/m c « 0.1 fm. 
Then it evolves into a J/^> whose mean size is quite larger, rj/^ rj 0.4 fm. 
Correspondingly, the absorption cross section, which scales as r 2 , increases 
by an order of magnitude. The expansion time is given by the uncertainty 
principle, 

Indeed, the produced cc dipole has a certain size, but no certain mass, and 
it takes time to resolve between the J/\l/ and the nearest radial excitation 
*£?'. There are in fact several time scales controlling the expansion process, 
the one given by Eq. (CQ) is the longest. 

A low energy dipole quickly expands to J/^f, while at high energy Lorentz 
time dilation freezes the initial small size for the time of propagation through 
the nucleus. So the nuclear medium becomes more transparent with rising 
energy, i.e. the effective break-up cross section decreases. 

To quantify this effect, let us consider a simplified equation based on the 
uncertainty relation describing the transverse expansion of a cc dipole moving 
with energy E Ec , 



dt E 5c 
Applying the uncertainty relation p T ~ l/r T , we get a solution 
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Here the initial separation squared of the cc dipole is fixed at the value 
(rj) ~ (5/rrjj, which deserves a discussion. At high energies the amplitude 



factorizes into the light-cone size distribution amplitude of cc fluctuations in 
a gluon, given by the modified Bessel function Kn(m c r T ), and the amplitude 
of cc dipole interaction with the target nucleon |l|. However, at low energies, 
when the time of charm production becomes as short as the proton radius, 
such factorization breaks down and the size distribution is poorly known. So 
in Eq. ([3]) 5 ~ 1, but not known more accurately (see more detailes in [l|). 
To evaluate the theoretical uncertainty we will try 5 = 0.5, 1.0, 2.0. 

Notice that the solution Eq. ([3]) is valid only for t < £/, when the quarks 
can be treated as free particles. At longer times the c-c interaction becomes 
important and affects the expansion process. A rigorous solution based on 
the path-integral technique is known [l|, |2J, but is more complicated. Here, 
for the sake of simplicity, we rely on the solution Eq. ([3]) assuming that 
t f > R Al i.e. E- cc > \R A {Ml - Mj /9 ). 

Due to color transparency |3( the dipole cross section in the small-r^ 
approximation has the form a a bs = Cr?, where the factor C(E 5c ) depends 
on dipole energy E 5c in the target rest frame. We can calculate the mean 
break-up cross section for a dipole of energy E 5c propagating and expanding 
along a path length L in a medium with a constant density, 

L 

a abs (L, E- cc ) = jjdl a aha {l) = C(E- CC ) f^- + -^ V (4) 

o 

We see that the mean break-up cross section is not a constant, as usually 
assumed, but rises with path length L and decreases with energy. The factor 
C(Ec C ) was calculated in J4J. For example, at xp = 0, and the energies of the 
experiments NA60 at CERN SPS and E866 at Fermilab, C(E Sc ) = 2.89 and 
3.14, respectively. 

Usually the nuclear ratio is evaluated with an oversimplified model as- 
suming that J/\l/ attenuates with a constant cross section <j a b s on the way 
out of the nucleus, 

R pA = -1— [ d% [I - e-^M] , (5) 

Aa abs J 

where a abs is treated as an unknown parameter fitted to data. As far, as 
we predicted the mean break-up cross section, Eq. (0J, we can calculate R pA 
and comparing with Eq. (JSj) extract a a b s - The result is plotted as function 
of energy in the left panel of Fig. [U This calculation was done with 5 = 
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Figure 1: Left: the break-up cross section fitted with expression ([5]) to the results of 
calculation with the absorption cross section Eq. (j4|) and S = 1 (solid curve), as function of 
energy. The bottom and top dashed curves show the theoretical uncertainty corresponding 
to variation of S = 1/2, 2 respectively Data are from fixed target experiments [5j,|6(. Right: 
ratio of the effective cross section Eq. ((4]) to the conventional mean one fitted with Eq. ([5]), 
as function of b at different energies of pA collision Ei a b — 158, 400, 800, 1200 GeV (from 
top to bottom). 



1 in Eq. 01]). The theoretical uncertainty is demonstrated by two dashed 
curves calculated with 5 = 1/2 (bottom) and 5 = 2 (top). Though with 
some uncertainty, our results explain well both the magnitude of a a i, s and its 
decreasing energy dependence [5|. 

Notice, that the extrapolation of a e ff up to the energies of RHIC, should 
be done with precautions. We remind that in this section we made an as- 
sumption about a short production time of a cc pair, which certainly breaks 
down at high energies. 

Since the effective break-up cross section rises with path length, it should 
be larger in central than in peripheral pA collisions. Indeed, in the right panel 
of Fig. [TJwe plotted the ratio of the L-dependent mean break-up cross section 
Eq. (@| to the one adjusted to the total J/^f cross section with Eq. (jSJ). The 
mean break-up cross section significantly exceeds the fitted one in central 
pA collisions, and underestimate it on the periphery. For this reason an 
extrapolation of nuclear effects from pA to AA with a constant break-up 



cross section cr a6s cannot be accurate. 

1.2. How long does it take to produce charm? 

Although the proper time of charm production is short, t* ~ l/2m c , in 
the rest frame of the nucleus, this time linearly rises with energy, 

tc - j^-. (6) 

Thus, if the energy of the produced J/^f is sufficiently high, E > 25 ( GeV) x 
L(fm), the effects of coherence become significant. This is a high twist 
shadowing in the process of a cc pair production by a projectile gluon. The 
cc is produced coherently in multiple interactions of the projectile gluon and 
the charm quarks with target nucleons. 

A good explanatory example is photoproduction of vector mesons on nu- 
clei. While at low energies the vector meson is photo-produced inside the 
nucleus and then attenuates through a half of the nuclear thickness, at high 
energies the vector meson appears as a Fock state of the incoming photon 
long prior the interaction and propagates through the whole nucleus. There- 
fore one expects a significantly stronger nuclear suppression at t c ^> Ra, than 
at t c <C R A [ll, 13]. Data for photoproduction of J/^f [8| and p mesons [9| 
nicely confirmed this prediction. 

Since the production amplitude is convoluted with the charmonium wave 
function, one can assume with a good accuracy an equal sharing of the total 
longitudinal momentum between c and c. Then, in the small-r^ approxima- 
tion the amplitude of cc production at the point with impact parameter b 
and longitudinal coordinate z inside the nucleus, averaged over the dipole 
size reads [Id], 



d 2 r T W 5c (r T )exp 



\c{E- cc )r 2 T (^T_(b,z)+T + (b,z] 



l+ l -C{E- cc )(rl) (^T_(b 1 z) + T + (b 1 z) 



(7) 



Here T.(b,z) = fl 00 dz f p A (b,z f )] T + {b,z) = T A (b) - T„(b,z), and T A (b) = 
T-(b, oo). Although the size distribution Wcc{tt) of produced dipoles has a 
complicated form, we assume for the sake of simplicity (more accurate calcu- 
lations will be published elsewhere) that it has a gaussian shape, W 5c (r T ) oc 



e -r T /(r T )^ with the mean value (rf.), which we estimated at (r^) = 6/m 2 . = 
0.1 fm 2 . Notice that due to color transparency the nuclear medium is more 
transparent than is expected in the Glauber model. Moreover, the ampli- 
tude Eq. (J7]) does not decrease exponentially with nuclear thickness, but as 
a power. 

Integrating the amplitude Eq. (J7J) squared over coordinates of the pro- 
duction point, one arrives at the nuclear ratio, which has the form, 

o = I/> fc IM ft* 

PA A J [1 + \C{E- CC ) (r|) T A (b)} [1 + l 2 C{E- cc ) (4) T A (b)} W 

At this point we can partially improve the small-r^ approximation in (JHJ) 
replacing C(E 5c ) (rf-) =>- °"$g( r T = ( r r))) where the dipole cross section has a 
saturated shape [2], and is somewhat smaller than C(E 5c ) r% at large ry. 

With Eq. (jSJ) we calculated the nuclear ratio Ra/ p {v) at i/i = 200 GeV, 
as function of J/^f rapidity y in the cm. of the collision, and its energy 

in the nuclear rest frame, E 5c = (y r s/2m N )./Mj, qi +p^e~ y . The results 

are depicted in Fig. [2] together with data from the PHENIX experiment 



11]. We see that the steep rise of the break-up cross section a 5c (r T} E 5c ) 
with energy (it triples from y = to y = 2) well explains the observed 
rapidity dependence of nuclear suppression. We should not continue our 
calculations far to negative rapidities, since the regime of long coherence 
length breaks down there. Besides, additional mechanisms, which cause a 
nuclear enhancement at negative rapidities, must be added. 

In the right panel of Fig. |2] we also present the impact parameter depen- 
dence of nuclear suppression for J/ty produced in proton-gold collisions at 
RHIC with different rapidities. As expected, the strongest dependence on 
rapidity comes from most central collisions. 

2. Gluon shadowing 

First of all, one should evaluate the kinematic condition for gluon shad- 
owing, t c ^ 9 > Ra, where tf 9 is the coherence time, or the lifetime of a ccg 
fluctuation in a gluon. This time can be related to the Ioffe time, as it was 



carefully calculated in [12] 



t ccg = J^ 5 /gN 

xm N 





Figure 2: Left: dashed curve presents nuclear suppression of J '/$ as function of rapidity in 
pA collisions calculated with Eq. jSJ. Solid curve is corrected for gluon shadowing. Data 
are for dAu collisions at y 7 ! = 200 GeV [ll| . Right: 6-dependence of the nuclear ratios for 
J/vP produced with rapidities y — 0, 1, 2, 3 in pAu collisions at ^/s = 200 GeV. 



0.1 was evaluated in 12 and found to be scale 



where the factor P g 

independent. Its smallness is caused by the large intrinsic transverse mo 



menta of gluons in hadrons, supported by numerous evidences in data [13 
Thus, shadowing for gluons onsets at smaller x < 0.01, than for quarks. 

In J/\l/ production at large X\ one should redefine x 2 =>- x 2 = x 2 /(l — X\) 
14| . Then the smallest x 2 = Mj^/sxii^l — x\) is reached at x\ = 1/2 
and equals to x 2 (mm) = 0.025 at at the energy of the E866 experiment, 
y/s = 40 GeV. This value is too large, so gluon shadowing has no contribution 
to the nuclear effects for J/\l/ production observed in the E866 experiment , 
as well as in other fixed target experiments 15 . 

Even at the energy y/s = 200 GeV, the values of x 2 are too large for gluon 

shadowing within the measured kinematics, (x 2 ) = e~ y J 2Mj ,^ 



(x 2 ) = e-y^2M^ + (p* T ) / 
where we use the cc invariant mass distribution predicted by the color singlet 



model 16|. With the measured (p\) = 4 GeV 2 the value of x 2 ranges from 
0.024 to 0.0033 within the measured rapidity interval < y < 2. We relied 
upon the results of the NLO analysis [l7| of DIS data, which suggest a ver 



weak gluon shadowing, in a good agreement with theoretical predictions [2|. 



The nuclear ratio corrected for gluons shadowing at Q 2 = 10 GeV 2 17 



is 



depicted in Fig. [2] by solid curve. We see that the effect of gluon shadowing 
is indeed vanishingly small. Even at the energy of LHC, y/s = 5.5 TeV and 
y = gluon shadowing according to 17J, |2| is extremely small, only 3% 
(x2 = 5.5 x 10 -3 ), and will be neglected in what follows. 

Notice that our explanation of the RHIC data is quite different from the 



description presented in [Uj. First of all, the charm quark shadowing was 
completely missed, and the naive formula ((SJ) with a fitted break-up cross 
section a abs was used. As we demonstrated, this formula is quite incorrect, 
especially for 6-dependence, even at low energies, where the approximation 
of short coherence length is reasonable. At the high energies of RHIC and 
LHC, such a formula is plain wrong. 

Further, the rapidity dependence of the nuclear ratio was prescribed 
in [ll| entirely to gluon shadowing, which had a rather large magnitude. 
Strangely, the authors referred to the same analysis |l7| as is used here, but 
they came up with a much stronger shadowing. Instead of the gluon shad- 
owing fitted in 
the authors of 



17] to data, they picked up a version called nDSg. However, 
171 ] warned that this version should not be used as a gluon 



PDF, since it "should be considered only as a mean to study variations on the 
gluon nuclear distribution" . For that purpose gluon shadowing was enforced 
to be large at small x contradicting data. 

A similar procedure was used in 18| . although for a different reason. The 
magnitude of gluon shadowing was fixed "by hands" at a large value at small 
x, otherwise the LO analysis in [18] would not have had any solution for gluon 
shadowing at small x. Because of this ad hoc input, the EKS shadowing is 
similar in magnitude, and is as reliable, as the nDSg. 

A strong gluon shadowing was reported recently in 



19 . Besides DIS 



data, this analysis includes data on hadron production in dA collisions at 
forward rapidities. Interpretation of this data is still controversial 15| , and 



this attempt to explain the observed nuclear effects entirely by coherence 
effects led to a gluon shadowing which significantly violates the unitarity 



bound 20 



3. Nontrivial transition from pA to AA 

At fist glance one might think that transition from nuclear effects in pA 
to AA collisions is straightforward: RaaQ), f) = R p a(t) x R p a{P — r), where 



nuclei collide with impact parameter b and J/^f is produced at impact param- 



eter r. Indeed, such a " data driven" procedure was used in [ll|, [21[ to predict 
the cold nuclear matter effects in nuclear collisions basing on measurements 
of 6-dependence of nuclear suppression in pA. 

3.1. Double- color- filtering 

The pA to AA transition, however, is not that simple. We illustrate 
this on the following example. If a cc dipole of transverse separation ry 
propagates through a slice of nuclear medium of thickness T A , its survival 
probability is S pA (tt) = exp(— Ct\Ta). Integrating over ry with the size 
distribution function W{j-t) oc exp[— -^/(rf,)] leads to (compare with ([7])), 

SpA = l + C(r* T )T A - (10) 

In the case of a central AA collision, according to the above recipe one should 
expect S AA (b) = Sp A . 

Let us, however repeat the above averaging over dipole size of S aa {tt) = 
Sp A {rx) = exp(— 2 x Cr^T A ). Factor 2 is here because the dipole attenuates 
simultaneously through both nuclei. The result of averaging over r? (left) 
should be compared with the conventional recipe (right), 

SAA = l + 2C(r*)T A * [1 + C(4)T A ] 2 (H) 

One can see that the two absorption factors are quite different, especially 
for C(r^,)T A > 1. The source of the difference is color filtering. Namely, 
the mean transverse size of a cc wave packet propagating through a nucleus 
is getting smaller, since large-size dipoles are filtered out (absorbed) with a 
larger probability. Such a dipole with a reduced mean size easier penetrates 
through the second colliding nucleus, compared to what would be in pA 
collision. The mutual color filtering makes both nuclei more transparent. 

Now we are in a position to perform realistic calculations for the nuclear 
suppression factor in AB collisions. Provided that the cc production occurs 
in the long coherence length regime for both nuclei, the nuclear suppression 
factor at impact parameter b reads, 

RAB(h) = 1 /,* KrYTMi-n 



x In 



Tab® J (A+-A A )(A+-A 

1 + A^ + A+)(1 + A+ + A B )- 



L(1 + A+ + A+)(1 + A A + A 



B) 



(12) 



where 
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C(3 



A(B)« 



T. 



A(B), 



(13) 
(14) 



-,A,B 



and E^ D are the energies of the cc in the rest frames of the nuclei A and B 
respectively. The result of calculation of Eq. (TT2|) is plotted by the upper solid 
curve in the left panel of Fig. |3J For comparison, the result of conventional 
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Figure 3: Le/t; Effects of double-color-filtering. J/^P suppression by the initial state 
interaction (ISI) effects in gold-gold collisions at %fs = 200 GeV as function of b. The 
upper and bottom pairs of curves (solid or dashed) correspond to y — and energies 
yfs = 200 GeV and 5.5 TeV respectively. Solid and dashed curves present the results 
at y = including and excluding the effect of double-color-filtcring, respectively. The 
dotted curve demonstrates rapidity dependence of the ISI effects at RHIC. It is calculated 
at y = 2 and is to be compared with the upper solid curve at y = 0. Right: Effects 
of boosted saturation scale. The upper and bottom dashed curves correspond to central 
gold-gold collisions at the energies *Js — 200 GeV, 5.5 TeV respectively. They demonstrate 
dependence on impact parameter r and are calculated in the same way as the solid curves 
in the left panel. Solid curves are calculated with the boosted saturation scale, which 
makes the nuclei more opaque for heavy dipolcs. 

calculations assuming simple multiplication of the suppression factors in the 
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two nuclei, is depicted by dashed curve. We see that the mutual color filtering 
makes the nuclei considerably more transparent. This effect should be more 
prominent for production of *&' and \. 

With Eq. f fl2|) we can trace the y-dependence of RaA. It turns out to 
be rather weak at the energy of RHIC, what obviously follows from the 
approximate linearity of y-dependence in pA depicted in the left panel of 
Fig. 121 However, at sufficiently large y, say y = 2, the condition of long 
coherence length breaks down in one of the nuclei. Then the cc dipole size 
is not frozen by Lorentz time delation, and the filtering in this particular 
nucleus is not effective any more. In this case the conventional multiplicative 
procedure is applicable, but the suppression factor in one nucleus (high E 5c ) 
should be calculated according to Eq. (|Sj), while in another nucleus (low 
E 5c ) one should do calculations for the short l c regime with the L-dependent 
absorption cross section Eq. 01]). The result of such calculation is plotted 
by the bottom solid curve in the right panel of Fig. [2j We see that the 
nuclear suppression at y = 2 is stronger than at y — 0. This happens due to 
disappearance of the double-color-filtering effect. 

3.2. Boosted saturation scale in AA collisions 

Another mechanism which breaks down the conventional multiplicative 
procedure for the transition from pA to AA is the mutual boosting of the 
saturation scale in AA collisions compared with pA [22]. It significantly in- 
creases the break-up cross section up to factor 1.5 and factor 3 at the energies 
of RHIC and LHC respectively. Correspondingly, the nuclear medium turns 
out to be much more opaque for cc dipoles in the case of nuclear collisions 



compared with the simplified multiplicative prescription of [11|, [21 . 

In the right panel of Fig. [3] we demonstrate the strength of this effect for 
central (b = 0) gold-gold collision as function of impact parameter r. The 
upper and bottom dashed curves corresponding to the energies of RHIC and 
LHC respectively, include the double-color-filtering effect, but exclude the 
saturation scale boosting, which is added to produce the solid curves. 

Thus, J/\l/ should be suppressed in AA collisions significantly stronger, 
than usually expected, and one should not interpret that as an anomalous 
suppression caused by final state interaction (FSI) with the dense medium. 
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4. Propagation of J/^f through a dense medium 

In the cm. of nuclear collision the nuclear disks passing through each 
other leave behind a cloud of radiated gluons creating a dense matter, which 
the J/\P propagates through. In this reference frame the J/^f full momentum 
is pt, which ranges from zero to several GeV in RHIC data. Such a low energy 
cc dipole develops the J/ty wave function pretty fast, during time tf < 0.5 fm 



16 , which is about the time scale of the medium creation. Thus, what is 



propagating through the medium is not a small cc dipole, but a fully formed 
J/ty. The mean dipole cross section is oj/^ = |C (rjm), where the factor C 



was introduced in fllj) and is known for a proton target from DIS data [23 



Its value for a hot medium is unknown, however, the factor C also controls 



broadening of a quark propagating through the medium [24j . So it is related 



to the transport coefficient q [25|, which is in-medium broadening per unit 
of length, C = q/2p. 

Therefore, the survival probability of J/^f produced at impact parameter 
f inside the medium has the form, 



R7a(?,pt) 



— exp 



rj/*> I dlq(f+l) 



(15) 



where \r + l\ 2 = r 2 + l 2 + 2rZcos0; Iq = vt ; and to — 0.5 fm. 

The transport coefficient depends on the medium density, which is func- 



tion of impact parameter and time. We rely on the conventional form [26 



{ J t n part (0,0)' l J 

where b and f are the impact parameter of the collision and of the point 
where the q is defined, respectively. The transport coefficient % corresponds 
to the maximal medium density produced at impact parameter r = in 
central gold-gold collision at the time t = t after the collision. We treat go 
as a adjusted parameter. 

The observed nuclear effects in J/ty production in AA collisions is inter- 
preted as a combination of FSI of J/^f in the dense medium Eq. ffT5| . and the 
initial state interaction (ISI) effects in production of J/^f caused by multiple 
interactions of the colliding nuclei. The latter was discussed above and in- 
cludes attenuation of the produced cc dipole propagating through both nuclei, 

12 



high twist shadowing of charm quarks, and leading twist gluon shadowing. 
In addition, gluon saturation in nuclei |4| leads to a considerable broadening 
of gluons, which causes a strong Cronin effect for J/ty. The details of cal- 
culations of the ISI effects can be found in 16 . The results for the nuclear 



effects in copper-copper and gold-gold are presented in Fig. HI The two plots 
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Figure 4: Nuclear ratio Raa for central copper-copper (full circles and squares, upper 
curve) and gold-gold (empty circles, bottom curve) collisions at y/s — 200 GeV as function 
of J /\& transverse momentum. The curves in the left and right p anels differ by calculation 
of the Cronin effect as described in text. Data are from 



differ by calculations of the Cronin effect. One can either make a shift in 
(p\) caused by broadening (left), or make a convolution of a nuclear- modified 
primordial transverse momentum distribution of the colliding partons with 
the known p^-distribution of J/ty (right). Both results are similar, except 
at large pt > 5 GeV, where no pp data are available, therefore no reliable 
prediction can be made. 

All effects are evaluated in a parameter free way, except the transport 
coefficient, which should be in the range of go ~ 0.3 — 0.5 GeV 2 / fm in order 
to reproduce data. This is close to the expected value % ~ 0.5 GeV 2 / fm 
25j, and more than order of magnitude less than was found from jet quench- 
ing data within the energy loss scenario 27]. Notice that the cold nuclear 



matter ISI suppression might have been underestimated in [16j. Our current 
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parameter- free estimate made in Sect. [3] results a significantly stronger ISI 
suppression, which is almost sufficient to explain RHIC data on J/ty produc- 
tion in central gold-gold collisions. This means that go may be even smaller, 
challenging the claim that a dense matter is created. 

5. Summary 

This talk highlighted several unusual features of J/ty production in pA 
and AA collisions, currently debated in the literature, which can be un- 
derstood taking a deeper look at the underlying dynamics. Since the wide 
spread interpretation of J/\l/ production off nuclei is grossly oversimplified, 
some improvements are proposed. 

• At the energies of fixed target experiments at SPS and Fermilab the 
break-up cross section for a cc dipole is subject to color transparency 
and is fluctuating during propagation through the nucleus. A simple 
model for the break-up cross section, which depends on energy and 
path length is developed. It well explains the energy dependence of 
the effective absorption cross section observed in data, and its absolute 
value. 

• At high energies of RHIC and LHC the charm production time becomes 
long, leading to a higher twist shadowing. This effect is of the same 
order as the attenuation caused by the cc break-up, and the magnitude 
of both is well fixed by DIS data from HERA. Data for J/ty suppression 
in dA collisions at RHIC are well explained without adjustment. 

• On the contrary, leading twist gluon shadowing is found to give no 
contribution to available RHIC data for J/\P suppression, and to be a 
rather small correction even at the energies of LHC. 

• Another effect, which makes transition from pA to AA model depen- 
dent is double color filtering. When a cc dipole propagates simulta- 
neously through the colliding nuclei, and one nucleus filters out large 
size dipoles, the reduced mean dipole size makes another nucleus more 
transparent. 

• Multiple interactions in the colliding nuclei lead to involvement of 
higher Fock states in the bound nucleons, which in turn enhance the 
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multiple interactions and bust the saturation scale further up. As a 
result, the nuclear medium becomes significantly more opaque for J '/\P 
in AA compared with pA collisions. 

• RHIC data for J/ty suppression in AA collision is well described com- 
bining the above effects and fitting the density of the produced hot 
medium, which is characterized by a transport coefficient q$. The found 
value of go is in good agreement with theoretical expectations, and is 
substantially smaller than what was extracted from jet quenching data. 
Thus, J '/\& production can serve as an efficient probe for the density of 
the created matter. 
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